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ABSTRACT 

We report the results of observations designed to investigate the spectroscopic signa- 
tures of dust clouds on the L2 brown dwarf Kelu-1. Time resolved medium resolution 
spectra show no significant evidence of variability in the dust sensitive TiO, CrH 
and FeH bandheads on the timescale of 1-24 hours. We do however report periodic 
variability in the psuedo-equivelent width of Ha consistent with the 1.8 hou r rota- 
tion period previously reported for this object ( [Clarke, Tinney fc Covey 2002| ). Near- 
contemporaneous /-band photometry shows evidence for non-periodic variability at 
the level of 2%. 



1 INTRODUCTION 



The brown dwarf Kelu-1 ( |Ruiz, Leggett fc Allard 1997| l 
is a L2 dwarf in the classification scheme of 
[Kirkpatrick et al. (1999)1 With an effective temperature of 
~1900K, refractory molecules such as TiO, CrH and FeH 
dominate the observed spectrum and play an important role 
in atmospheric physics. [Clarke, Tinney fc Covey (2002) | 
(hereinafter CTC) have discovered periodic variability 
in photometric observations centred on the complex of 
molecular bandheads at ~8600A. They provide several 
possible explanations for the variability, including: 1) dust 
cloud imhomogeneities modulating the surface brightness 
as Kelu-1 rotates (with an implied period of 1.8 hours), 
2) A close substellar binary inducing ellipsoidal variability 
(implying an orbital period of 3.6 hours). In this paper, we 
report the results of a search for spectroscopic variability of 
Kelu-1 which answers this question. 

Several spectroscopic searches have previously 
been made for variability in L and T dwarfs. 
[Nakajima et al. (2000)[ found evidence for possible variabil- 
ity in the near IR water lines o f the T dwarf SDSS1624-f 00, 
and [Kirkpatrick et al. (2001)"] detected changes in the 
~8700A CrH feature of the L8 dwarf Gl 584C. W ith 
very similar aims to this study, [Bailer- Jones (2002)[ has 
carried out near IR spectroscopic monitoring of the variable 
LI. 5 brown dwarf 2MASSW J1145572+231730 (2M1145). 



* Based on observations obtained at the European Southern Ob- 
servatory, Paranal, Chile (ESO Programme 68.C-0381) 
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Photometric observations of 2M1145 show variability, but 
no evidence for periodicities, which Bailer- Jones fc Mundt 
(2001a,b) claim as evidence for surface features evolving on 
the timescale of the rotation period. Spectroscopic obser- 
vations spanning 54 hours do not show any evidence for 
variability. Bailer- Jones places an upper limit of 10-15% on 
the covering fraction of clear holes in a dusty photosphere 
(similar to our model 3 in EIJ. 

In this paper we present the results of combined op- 
tical photometric and spectroscopic observations of Kelu- 
1 designed to test the causes of variability proposed by 
[Clarke, Tinney fc Covey (2002)[ Section 2 describes the 
data acquisition and reductions. Analysis of the resulting 
spectra and lightcurves are presented in section 3, and in 
section 4 we develop toy models to investigate our observa- 
tions. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Spectroscopy 

Kelu-1 was observed on two consecutive nights (2002 Febru- 
ary 13 and 14 UT) with the F0RS2 instrument on VLT 
UT4 (Yepun). Both nights were photometric, with subarc- 
second seeing. A log of observations is given in table Q The 
3001 grism and l"slit were used, providing a spectral res- 
olution of ~13A over the range 6300-11500A. The OG570 
order blocking filter was also used to remove second order 
light. The slit was aligned such that Kelu-1 and a brighter 
comparison star 25 arcsec to the south-west could be ob- 
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Date 


UT time 


Seeing 


Sky conditions 


2002 Feb 13 


05:10-07:20 


0.75-0.9" 


Photometric 


2002 Feb 14 


05:22-07:15 


0.75-1.0" 


Photometric 



Table 1. Log of FORS2 spectroscopic observations carried out 
on VLT UT4. 

served simultaneously. Figure Q shows the alignment of the 
slit, and identifies Kelu-1 and the comparison star. The two 
background stars seen near the slit in figure^are several arc- 
seconds from the slit, and there is no detectable scattered 
light from them in our FORS2 spectra. 

Each night's observations consisted of a two hour se- 
quence of 480s integrations (for a total of 13 spectra per 
night). This gave us complete phase coverage of the 1.8 hour 
rotational period suggested by CTC, and provided enough 
signal to noise (~70-80 per pixel; ~160-180 per resolution 
element) to detect small changes between consecutive spec- 
tra. In addition to the science exposures, the ESO calibra- 
tion plan provided lamp flats and observations of the spec- 
trophotometric standard CD-32 9927 Pamuy et al. 1994| |. 
Observations of the standard were taken about 1 hour after 
the science observations on night 2. HeAr arclamp exposures 
for wavelength calibration were also taken at the end of the 
2nd night. 

Observations were reduced with IRAF^. After subtract- 
ing a bias frame, and correcting pixel-pixel variations and 
fringing with a normalized lamp fiat, the spectra were ex- 
tracted with APALL task using the optimal extraction 
method of |Horne (1986)| The spectra were wavelength cal- 
ibrated via the HeAr arc spectra. To keep the spectra as 
"raw" as possible, we did not interpolate them onto a linear 
scale. The night sky lines show that the wavelength solution 
is stable during each night, but changes by ~4A between the 
nights. The spectrophotometric standard was then used to 
correct the illumination function of the instrument. Figure|5| 
shows the averaged spectrum of Kelu-1 from the 1st night. 

2.2 Photometry 

7-band photometry of Kelu-1 was obtained in service mode 
with the EMMI instrument at the NTT on the nights of 
2002 Feburary 17 and 19 UT. ESO fiher #610 was used, 
which is slightly redder than the standard Cousins Ic filter. 
For the rest of this paper, we have used the ESO #610 fil- 
ter profile when discussing the "/" filter. The CCD chip was 
windowed to give a 3.2' x3. 2' field of view. To further reduce 
readout overheads we employed 2x2 pixel binning, giving a 
pixel scale of 0.54"/pixel. The observations consisted of con- 
secutive 120 second exposures over two hours, with around 
forty frames obtained each night. 

Data were reduced in the standard fashion with IRAF, 
subtracting a bias frame and dividing by a normalized 
dome fiat. Aperture photometry was then performed on 
the target and a selection of comparison stars with the 
APPHOT package within IRAF. Diflterential hghtcurves 

^ The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by NOAO, which is operated by AURA, Inc., under con- 
tract to the NSF. 
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Figure 1. Slit position on the sky of the FORS observations 
showing the position of Kelu-1 and the comparison star. 
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Figure 2. Spectrum of Kelu-1 with important atmospheric fea- 
tures denoted. The insert panel shows a close up of the molecular 
complex at ~8600A 

(target— comparison) were then constructed in the man- 
ner described by [Clar ke, Oppcnhci mer fc Tinney (2002)| 
Briefiy, we produce an mean comparison star from an ensem- 
ble of non-variable stars in the field, and it is this "mean" 
that we refer to as the comparison star in H3.ll 

3 RESULTS 

3.1 J-band photometry 

Figure 13 shows the differential 7-band lightcurve of Kelu-1 
minus comparison star (upper panel) and comparison minus 
check star (lower panel). Table |5| gives the statistics of the 
hghtcurves. The rise of 0.02 mag in the Kelu-1 lightcurve 
from 4.34 to 4.36 (~30 minutes) is significant at the 4-sigma 
level. The comparison minus check lightcurve remains con- 
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Lightcurvc 


Mean 


Standard Dev 


Average 


Dev 




night 1 


night 2 


night 1 


night 2 


night 1 


night 2 


Kclu-1 — Comp 


-0.0082 


0.0117 


0.0069 


0.0054 


0.0049 


0.0042 


Check — Comp 


-0.0013 


0.0012 


0.0078 


0.0065 


0.0057 


0.0045 



Table 2. Statistics for Kelu-1 NTT differential photometry. The mean for each night is mean magnitude of that night's data relative to 
the mean magnitude of the whole dataset. Standard deviation is defined as ^^J (xj — x)^ and the average deviation is defined as 



slant to 0.002 mag during this time. This rise represents a 
dimming of Kelu-1. 

Periodogram analysis of the data does not show evi- 
dence for any significant periods, but the analysis is heavily 
compromised by the window function of the data. They are 
however, inconsistent with the 1.8 hour (and 3.6 hour) pe- 
riod found in the 2000 photometry. As discussed later ( t|3.3|| . 
there is evidence from the spectroscopy to support the in- 
terpretation of the 1.8 hour period as the rotation period. 
The 7-band and R6 fiher used by CTC (A^en = 8580A, 
AA = 410A) should react similarly to inhomogeneities in 
the atmosphere (see Figure [TTI and H4.1II . so we must as- 
sume that the atmosphere of Kelu-1 has evolved signifi- 
cantly between March 2000 and February 2002. It may be 
that in March 2000, Kelu-1 had a significant long lived 
(>25 rotations) surface feature, which has since dissipated. 
In February 2002, if stable features were present, they 
were smaller and affected the lightcurve at less than the 
0.5% level. The 0.02 mag dimming may be due to more 
rapid evolution of a surface feature, as has been claimed 
for other brown dwarfs lIBailer- Jones fc Mundt 2001al 
lEall er-Jones fc Mundt 2001 bll . 



3.2 Dust senstitive spectroscopic features 

At the temperature of Kelu-l's photosphere (Tgjj~1900K) 
gas phase molecules begin to condense into solids (i.e. dust). 
The removal of species like Ti, V, Fe and Cr from the gas 
phase into solids results in the weakening of TiO, VO, FeH 
and CrH bands. If Kelu-l's photosphere has significant spa- 
tial dust imhomogeneities, we may expect to see the signa- 
ture in these lines (see tl4.1ll . The complex of lines around 
8700A (TiO, CrH and FeH) is particularly interesting, as 
these lines dominate the photometric band in which CTC 
detected variability. 

Figure|3]shows the time series spectra of the 8200-8900A 
region, and does not show any obvious evidence for variabil- 
ity. To make any small differences more apparent to the eye, 
we have plotted the difference specta of this region in figure|S] 
An average spectrum has been constructed for each night, 
and this has been subtracted from each individual spectrum. 
In both plots, each spectrum has been offset for clarity, with 
time increasing from the top of the diagram. The offsets be- 
tween consecutive spectra are constant however, and do not 
represent the actual time between observations. There does 
not appear to be any evidence for variability of the molecular 
bands either between or during each night. 

To test for subtle changes in the absorption depth, we 
have measured the flux ratio between narrow bands stradling 
the molecular bandheads. These bands are shown in right- 



Band 


Acen 


AA 


TiOl 


8395 


50 


Ti02 


8485 


50 


CrHl 


8580 


40 


CrH2 


8628 


16 


FeHl 


8670 


20 


FeH2 


8710 


20 



Table 3. Bands used to measure the depth of molecular absorp- 
tion features. 

hand panel of figure 01 and described in table |21 The flux 
through each band was summed, and the molecular band 
index is defined as the ratio of the fluxes, i.e. TiO index 
= Ti01/Ti02. To check for instrumental or atmospheric 
effects, we also measured band indices for the comparison 
star, although no molecular absorption is seen in this ob- 
ject. Figure |n] shows the TiO, CrH and FeH indicies during 
the observations for Kelu-1 (upper panels) and the compar- 
ison star (lower panels) . Errors are carried through from the 
error estimates produced by APALL. The band indices for 
all these molecules are constant to within ±2%. The most 
significant indication of variability is in CrH between the 
nights. This is significant at the 3(t level, but the fact that 
similar effects are seen in the comparison star leads us to be- 
lieve we are not observing effects intrinsic to Kelu-1. The la 
upper limits on variability (peak-to-peak) of the molecular 
band indices are; TiO<1.4%, CrH<4%, FeH<2.2%. 

3.3 Hq Variability 

The Hq data show evidence for significant variability on 
even a cursory examination. Figure |7| shows a close up the 
Ha line from all the spectra. Each spectrum has been offset 
for clarity, in the same fashion as figure 2] 

We have measured the psuedo equivelent width (PEW) 
of the Hq line in all the spectra with the splot task in IRAF. 
The results are plotted in figure |S] as PEW(Hq) as a func- 
tion of time. It is clear from this plot that the PEW(Hq) 
changes with time. Periodogram analysis reveals a period in 
the range 1.5-2.5 hours - consistent with the 1.8 hour pe- 
riod photometric period reported by CTC. The broad peak 
in the periodogram is due to the limited time coverage on 
each night. Figure|U]shows the PEW(Hq) folded on a period 
of 1.8 hours. The central wavelength of the Hq line does not 
change significantly throughout the observations. 

The mean level of PEW(Hq) is ~3.0A (in emis- 
sion), which is slightly larger than, although com- 
parable to previously measured values (e.g. 1.5-2.0 A; 
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Figure 3. /-band lightcurve of Kelu-1 on 2002 February 17 and 19 UT. The upper panel shows Kelu-1 minus comparison star, and the 
lower panel shows comparison minus check star. Each lightcurve has been mean subtracted (using data from both nights). The thick 
bars show the photometry averaged into 10 minute bins. 
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8200 8400 8600 8800 8200 8400 8600 8800 



Wavelength (A) Wavelength (A) 

Figure 4. Time series spectra of the molecular lines of TiO, CrH and FeH at 8400— 88OOA. Time increases downwards, but the plotted 
gap docs not represent the actual time between observations. In the right panel (night 2), we have marked the bands used to measure 
the strength of the absorption features (Table 1^ 
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Figure 5. Difference spectra of tlie molecular lines of TiO, CrH and FeH at 8400— i 
night's spectra (upper most spectrum) subtracted. The increased noise redward of 
and lower S:N in this region (~75% of the SN at 86OOA). 



Each spectrum has had the average of all the 
is due to a combination of telluric variability 



IBasri 200HKirkpatrick et al. 1999| l. We have also mea- 
sured the flux in the Ha hne with SPLOT. The mean 
Ha hne flux during our observations was 0.8 ± 0.4 x 
10~^® erg/s/cm^, comparable with the value reported by 
[Ruiz, Leggett fc Allard (1997)| To check the flux calibration 
of our spectroscopy, we have calculated the flux through the 
I bandpass. This is ~75% higher than expected (given Kelu- 
I's / band magnitude of 16.8, [Ruiz, Leggett fc Allard 1997| l, 
but the discrepancy is probably due to overestimate of flux 
towards the red end of the spectrum, where the standard 
star is very faint. We therefore assume a conservative error 
of 50% in our line flux (included above). 



3.4 Radial velocity changes 

To look for changes in Kelu-1 's radial velocity, we have mea- 
sured the wavelength of the strong neutral atomic lines of 
Rbl (7800, 7948A) and Csl (8521 A). These lines lie in rel- 
atively "clean" areas of the spectrum and, unlike the close 
doublets of KI and NaT, allow simple single proflle fitting. 
We used the splot task in IRAF to do this. To check for 
instrumental drifts, we measured the central wavelength of 



several lines in the comparison star. For each line, we con- 
vert the measured wavelength to a velocity relative to the 
mean wavelength of the line in all spectra. All the lines are 
then averaged to give a radial velocity for each spectrum. 
Figure [THl shows the radial velocity curves for Kelu-1 (up- 
per panel) and the comparison star (lower panel). The same 
trend is seen in both stars, indicating the observed velocity 
changes are instrumental effects. This is most likely caused 
by drifts in the centering of the stars on the slit, which was 
typically wider than the FWHM. The scatter on the velocity 
measurements of Kelu-1 is ±10kms~^. This provides an up- 
per limit on radial velocity variations, and rules out a close 
companion with Msini >10Mj,jp^. We can therefore reject 
the binary hypothesis for Kelu-l's variability proposed by 
CTC. Note also that the lack of a 1.8 hour period in the I- 
band photometry also rules out a companion as the source 
of variability, as a companion would produce consistent and 
repeatable variability. 



Assuming a mass of O.O65M0 for Kelu-1 
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Figure 6. Changes in tlie molecular band strengths during the 
observations for Kelu-1 (upper panels) and the comparison star 
(lower panels). The definition of band strength of given in the 
text. 
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Figure 7. Close ups of the Ha line. The first night is shown on 
the left, and the second night on the right. The spectra have been 
offset for clarity, with the earliest spectra in each night top of 
each panel. 



4 DISCUSSION 

4.1 Limits on the inhomogeniety of Kelu-l's 
atmosphere 



[Allard et al. (2001)| provide theoretical spectra for models 
relating to the two limiting cases of dust formation. The 
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Figure 8. PEW(Ho) during the FORS observations. The upper 
panel shows night 1, and the lower panel night 2. Typical errors 
on the PEW(Ha) are O.SA, as shown by the vertical line in upper 
panel. 
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Figure 9. PEW(Ha) folded on the period of 1.8 hours reported 
by CTC. This period is consistent with a periodogram analysis 
of the PEW(Ha) data itself. 
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Figure 10. Radial velocity measurements of Kelu-1 (empty cir- 
cles) and the comparison star (filled squares). The times of the 
comparison star points have been slightly offset for clarity. Each 
meaurement is the average of several lines, and the error bars are 
given by the standard deviation. 
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DUSTY model assumes that dust forms and remains sus- 
pended evenly throughout the photosphere, analogous to 
small grains supported in thermal updrafts. Alternatively, 
in the COND model, dust particles form but immediately 
settle below the photosphere - a good model for grains which 
are too large to be supported by convective updrafts. In re- 
ality, the detailed physics of dust formation in substellar 
atmospheres is much more complex, being a coupled system 
involving chemisty, atmospheric dynamics, radiative trans- 
fer and cloud formation. In particular, gravitational set- 
tling of dust grains and thermal structure in the atmosphere 
will lead to a vertically stratified distribution of dust. Sev- 
eral groups fAckcrman & Mar le y 200 1] [Cooper et al. 2002| 
[Tsuji 2()02nMarley et al. 2002J have made more detailed at- 
tempts to model dust formation in substellar atmospheres. 
All models to date, however, are static, and do not treat any 
horizontal inhomogeneities in atmospheric structure. 

To make a crude estimate of the effects a "patchy" 
photosphere would have on the emergent spectrum, we can 
combine the DUSTY and COND models. The overall spec- 
trum of Kelu-1 is best fit by a DUSTY model at 1900K 
(Baraffc ct al . 199811 . so we assume that inhomogeneities are 
best modeled by COND "holes" in an otherwise DUSTY 
photosphere. By definition, holes will allow us to see into 
deeper, hotter regions of the photosphere. With this in mind, 
we have constructed several models; 

Model-1 2100K COND holes representing patches of ef- 
fcient dust settling which allow us to see deeper into the 
atmosphere. 

Model-2 2100K DUSTY holes representing patches 
which allow us to see a deeper level of the atmosphere, 
but where dust is still significant. 

Model-3 1900K COND holes representing patches of ef- 
ficent dust settling at the same temperature as the rest 
of the photosphere. 

In all cases, the 1900K DUSTY spectrum represents the ma- 
jority of the photosphere. To model the effects of imhomo- 
geneties in the atmosphere, we construct a linear combina- 
tion of the two spectra as; 

/(^)mix = (1-0 - ^^)./(A)phot + Ff Whole (1) 

where /(A)pjjQ^ and /(A)jjQjg are the model spectra repre- 
senting the normal and "hole" regions of the atmosphere, 
and F is the covering fraction of holes. We therefore gener- 
ate a "mixed" spectrum, for which we can measure molec- 
ular band indices in exactly the same way as for the real 
data ( ^n\ . We also measure the photometric bands I 
(ESO #610) and R6 (as used by CTC; Acen = 8580A, 
AA = 410A). 

Figure fTTI shows the efi'ects of changing the "hole" cov- 
ering fraction for all three models. The band indices are 
normalised to the purely dusty case (no holes, F = 0), so 
the points give the change in band indices we would observe 
if the photosphere changed from having no holes, to having 
a given covering fraction. The effect of small changes in F 
is roughly linear, so the plots are also valid for a changes 
between two non-zero covering fractions. 

The horizontal shaded area shows the R6 photometric 
amplitude reported by CTC (1.2±0.1%). The most striking 
feature of these plots is the incredibly small changes in cov- 



ering fraction (vertical shaded area) required to reproduce 
the observed photometric variability; less than 2% "holes" 
in all three models. Comparing with the molecular band in- 
dices, our upper limits on variability (TiO<1.4%; CrH<4%; 
FeH<2.2%) are entirely consistent with these covering frac- 
tions for all models. It should be noted however, that these 
limits represent the difference in covering fraction between 
one side of the brown dwarf and the other, the atmosphere 
may be far more inhomogenous on small scales. 

Comparing the models reveals interesting differences 
between them. The CrH index appears to be much more 
sensitive to dust properties (models 1 & 3) rather than 
temperature (model 2). The TiO index offers a potential de- 
scriminant between models 1 & 3, becoming stronger with 
covering fraction in model 1, and weaker in model 3. None of 
the molecular bands are strongly effected in model 2, which 
depends purely on temperature. The fact that the scatter 
in our measurement of CrH is larger than TiO and FeH is 
more in line with models 1 and 3 than model 2. However, the 
most probable explanation of differences between the mod- 
els is that our toy models are too primitive to throw any 
light onto these subtle effects. It is likely a full 3-dimensional 
model of a brown dwarfs atmosphere will be required to in- 
terpret the results of future spectroscopic studies. 

4.2 The chromosphere of Kelu-1 

The Hq emission we see is probably caused by ionised hy- 
drogen gas in Kelu-1 's chromosphere, analogous to other 
low mass stars (e.g. jGizis ct al. 2000J. The energy re- 
quired to heat the gas may come from magnetic heat- 
ing ( |Mohanty et al. 20'02| |, or wave heating llYelle 200011 . al- 
though calculations of the energy released from both are 
problematic. Combined with magnetic fields above the sur- 
face of Kelu-1, this could lead to inhomogenous patches of 
hot gas above Kelu-1 and produce the variable Ha emis- 
sion we observe. Similar periodic Ha variability has been ob- 
served from other stars. [Fternandez fc Miranda (1998)| show 
the variable Ha emission from several weaklined T-Tauri 
stars is correlated with photometric variability, which they 
interpret as chromospheric magnetic loops corotating above 
cool magnetic spots in the photosphere. Similar anchored 
magnetic loops could exist above Kelu-1, leading to Hq vari- 
ability with the same (rotation) period as photospheric vari- 
ability. 

The presence of magnetic fields and a chromo- 
sphere do not necessarily contradict previous arguments 
against magnetic spots in the photospheres of ultra cool 
dwarfs fB ailcr- Jones & Mundt 2001a Gclino ct al. 2003 
|Mohant~et al. 2002.. Bailer- Jones 2002.1 . Densities in the re- 
gion of the photosphere are much higher, and the atmo- 
sphere there is much more neutral. This corresponds to 
very low magnetic Reynolds numbers (Gclino ct al. 2002J, 
or high electrical resistivities jMohanty et al. 2002} . In the 
region of the photosphere, the magnetic field is decoupled 
from the atmospheric fluid. A magnetic field may be gener- 
ated deeper in the atmosphere, where the temperature and 
ionisation fraction rise, and the magnetic field is well coupled 
to the matter. 

The Hq line flux we measured in il3.3l corresponds to 
an Ha luminosity of 3.7 ± 1.8 x 10^* erg/s/cm^ at Kelu-l's 
distance of 19.6 pc lIDahn et al. 20021 . In the standard mea- 
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Figure 11. Changes in band strength with covering fraction for all three models. Model-l upper left, Model-2 upper right and Model- 
3 lower left. The shaded areas in each plot represent the amplitude of the R,6 variability discovered by CTC (1.2±0.1%), and the 
corresponding surface covering fraction. A covering fraction of represents a photosphere with no "holes" . 



sure of chromospheric activity, Kelu-1 has log(L(Ha)/Lj,Qi) 
of —5.35 ± 0.5, where our error estimates are quite conser- 
vative ( H3.3t . This value is in line with what we would ex- 
pect for a L2 dwarf from the work of |Gizis et al. (2000)| and 
[Mohanty et al. (2002)] 



5 CONCLUSIONS 

We have presented here the first high signal-to-noise phase 
resolved optical spectroscopy of a brown dwarf with a known 
rotation period. We have detected a rotational modulation 
of the Ha line, consistent with the 1.8 hour rotation period 
reported by |Clarke, Tinney fc Covey (2002)] Radial velocity 
measurements reject the possibility that photometric vari- 
ability is caused by a close companion to Kelu-1. Photom- 
etry shows that the atmosphere of Kelu-1 has significantly 
evolved between March 2000 and February 2002, when no 
periodic signal can be detected at the level of 0.5%. Kelu- 
1 does however show evidence for a 2% dimming over ~30 
minutes. We have also placed upper limits on the scale of 
possible surface features inducing inhomogenous dust for- 
mation across the surface, and made a primitive attempt to 
investigate the effects of different atmosphere models. On 
the balance of observations, it seems that Kelu-1 is a single 
brown dwarf, with a typical level of chromospheric activ- 



ity, and a relatively homogenous, although evolving, atmo- 
sphere. 
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